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We provide a path integral treatment of the decoherence process induced by a heat bath on a 
single particle whose dynamics is classically chaotic and show that the decoherence rate is given by 
the Lyapunov exponent. The loss of coherence is charaterized by the purity, which is calculated 
semiclassically within diagonal approximation, when the particle initial state is a single gaussian 
wave packet. The calculation is performed for weak dissipation and in the high temperature limit. 
This situation allows us to simplify the heat bath description to a single random potential. Although 
the dissipative term is neglected in such approach, the fluctuating one can be treated phenomeno- 
logically to fit with the above regime. Our results are therefore valid for times shorter than the 
inverse of the dissipation rate. 
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I. INTRODUCTION 

A well-established way of making classical features to 
appear in quantum systems is to study them when the 
actions involved are much greater than Planck's con- 
stant. This semiclassical approach has allowed us to ex- 
plore the quantum-classical transition in a very interest- 
ing class of systems namely those which are classically 
chaotic. In different contexts, several results have shown 
how the classically chaotic dynamics has an influence on 
the quantum behavior However, it has been shown 
experimentally that strictly quantum features such as 
superpositions between macroscopically different states 
are still present even when actions are bigger than Ti. 
In Ref.[2j, interference patterns measured in Josephson 
junctions show the presence of superpositions between 
current states, each of those corresponding to the mo- 
tion of millions of Cooper pairs. Therefore, besides the 
semiclassical regime, a decoherence mechanism must be 
responsible for the emergence of pure classical behavior 
where such superpositions are absent [3|. 

The decoherence process is one of the effects induced 
by the coupling to an environment and the path integral 
description of linear systems interacting with a heat bath 
is very well established 043- Nevertheless, for nonlinear 
systems there is no general path integral approach to such 
problem and specific kinds of approximations have been 
developed @. At first glance, it can be argued that de- 
coherence happens so fast that does not matter whether 
the dynamics is linear or not Q. On the other hand, 
calculating a decoherence rate in such regime contains 
the assumption that the isolated dynamics of the sys- 
tem of interest is not relevant at all. That is true only 
if the coupling between system and heat bath is weak 
and this is not the regime we are interested in here. 
Once the nonlinearity is relevant, different semiclassical 
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approaches are available [9l-lllj either within the path in- 
tegral or the master equation framework. 

In this context, a very interesting result was conjec- 
tured for the first time in Ref.[l2|]. There, the authors 
proposed theoretically that a classically chaotic system 
could lose coherence at a rate given by the Lyapunov 
exponent, i.e. a rate completely independent of the envi- 
ronment parameters and related to the classically chaotic 
dynamics. A series of works appeared afterwards trying 
to implement that proposal in solid basis (focus- 
ing also on the quantum-classical correspondence aspect 
of the problem). A first approach was developed based 
on the resolution of a master equation for the Wigner 
function of the classically chaotic system. It was shown 
analytically [lij and numerically [l6[ that the entropy 
production rate obtained from that approach is propor- 
tional to the Lyapunov exponent. The analytical calcu- 
lation was however obtained under the following approx- 
imations: the chaotic dynamics is taken into account up 
to its linear regime, the friction from the heat bath and 
the quantum corrections to the classical evolution of the 
Wigner function are neglected. Despite of this, such re- 
sults were suported by more rigourous calculations con- 
sidering a hyperbolic linear system [l8j . 

A second approach was developed using path integral 
methods and it is known in the literature by Loschmidt 
echo (LE) [I9I - I21I ]. The LE is the probability to recover 
a certain initial state when the Hamiltonian used in the 
forward time evolution is slightly perturbed for the time 
reversed propagation of the final state. Considered as 
an indirect measure of decoherence [H|-[25j], the LE was 
shown to decay with the Lyapunov exponent when the 
system is classically chaotic. On the other hand, the LE 
deals with isolated systems only and therefore it has been 
understood as a stability measure of them, as proposed 
originally by Peres 26]. Thus, a path integral formu- 
lation of a classically chaotic system interacting with a 
heat bath still is a very challenging problem. 

The main goal of the present work is to provide a path 
integral derivation of the result conjectured in [l2| avoid- 



2 



ing the LE problems mentioned above. In Sec. II the 
treatment of the heat bath is described and the regime 
considered is presented. In Sec. Ill we consider the purity 
as a measure of decoherence and it is calculated through 
semiclassical methods for arbitrary initial states. In Sec. 
IV it is shown the role played by correlations between 
pairs of classical trajectories in the decay of the purity 
and how they give rise to a decoherence rate given by the 
Lyapunov exponent when the initial state is a gaussian 
wave packet of arbitrary width. Finally, conclusions are 
presented in Sec. V. In the appendix we present a brief 
comparision between LE and purity calculations pointing 
out their differences. 



II. DESCRIPTION OF THE HEAT BATH 

We consider here a single particle whose dynamics is 
classically chaotic as, for example, electrons in a clean 
ballistic conductor [27]. Our goal is to describe its dy- 
namics when it is coupled to a heat bath, focusing on 
the decoherence process. We start from the path integral 
formulation in terms of Feynman- Vernon influence func- 
tionals [28| with a Lagrangian given by L = Ls+Lj+Lb- 
The Ls term describes the isolated dynamics of the par- 
ticle which we assume to be a classically chaotic billiard- 
type of confinement. The terms L\ and Lb, interaction 
and heat bath Lagrangians respectively, are assumed to 
be given by Caldeira-Leggett model [IJ with a linear cou- 
pling between particle and bath coordinates. For initial 
states of the form p = ps ® pb, where ps and ps are 
the initial density matrices of the particle and the bath 
respectively, the time evolution of the reduced density 
matrix ps is given by 28] 

Ps (x/ , y/ ; t ) = J dXi dy l J(x/ , y f , x, , y t ; t ) p s (x 4 , y, ; ; 0) , 

(1) 

where the propagator J, obtained after tracing out the 
heat bath degrees of freedom with ps chosen as a thermal 
state, is given by the following path integral 

J(x /! y / ,x i ,y i ;t)= [vxVyc^ff^, (2) 



For an Ohmic spectral densityand in the high temper- 
ature limit, S e ff is given by [34a] 

Se//[x,y] = 

(3) 

<S [x,y] +?7i7$ d [x,y] +^ $/[x,y], 

S [x,y] = 5 [x]-5 [y], S [q\ = J* ds L s [q.(s), q(s)] gives 
the isolated particle dynamics and the other two terms, 
$d and $/, lead to dissipation and decoherence respec- 
tively. They are local quadratic functionals of the paths 
x and y without any other prefactors involving the par- 
ticle mass 771, the friction constant 7, h or the product 



fcgT of the Boltzmann constant and the heat bath tem- 
perature 043 ■ I n a regime of very weak dissipation and 
sufficiently high temperature, the term $d is negligible 
at least for time scales much shorter than 1/7. In this 
situation, the effective influence of the heat bath can be 
modeled in a much simpler way. 

Let us consider for a moment the following Lagrangian 
L = Ls — q • f(i) where q is the position vector of the 
system Ls and f(t) is a stochastic force. One can also 
construct the time evolution of the reduced density ma- 
trix ps in this case but it would depend on f(t). Thus, 
assuming that f(t) is gaussian distributed with an av- 
erage value (f(t)) = 0, we can calculate the following 
averaged propagator [29| 

J(x f ,y f ,x i ,y i ;t) = / VtP[f\ [ VxVy e*^ f ],(4) 



where P[f] is a gaussian functional and 

S[x, y; f] - S [x, y] - J ds f (a) ■ (x(s) - y(,s)) . (5) 

The average over f(t) leads to the effective action 
S e ff = 5 G [x,y] + iS f [x,y] in which Q 



<5/[x,y] = 



ds / ds 
Jo 



'(x( S )-y( S ))M( S - S ')(x( S ')-y(s')) 

(6) 



and M(s — s') is 2 x 2 matrix (since we are consider- 
ing a two dimensional system) whose elements are the 
correlation functions of the stochastic force 



Mi, j (s- S ') = i(/ i ( S )/ i ( S ')), 



(7) 



with /j(t), i = 1,2, being the f(t) components. Since we 
are free to choose the spectra of (fi(s)fj(T)), we take it 
as 

^(Ms)fiW)) = ^^5^8-8'). (8) 

where 5ij is the Kronecker delta. This choice leads to 

S ef f [x, y] = S [x, y] + » 4m ^ T $/ [ X; y ] . (9) 

Eq.© is identical to ([3]) except for the dissipative term. 
Therefore, we can model the influence of the heat bath in 
the regime of very weak dissipation and high temperature 
just by the fluctuating time dependent potential shown 
in ([5]) and by replacing the propagator ([2]) by the one in 
©■ 



III. SEMICLASSICAL PURITY 

The next step then is to write a semiclassical expres- 
sion for Q, which can be written in terms of Feynman's 
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propagators 



•7(x/,y/,Xi,yi;i) = 

J VfP[f]K( Xf , Xi ,f;t)K*(y f ,y h f;t) (10) 



where if (x/,X£ } f;i) = / 2? X e 4S [ x ' f l/' i , if* is its complex 
conjugate and 5[x, f] = 5 D [x] — f Q dst(s) ■ x(s). 

We assume that the coupling to f is classically small in 
order that only the actions, i.e. the phases, are affected 
while the trajectories given by Ls remain unchanged. 
Under this assumption, we replace K by the semiclassical 



Van Vleck formula [l] for two dimensional systems 
K sc (xf,Xi,f;t) = 

(dfi) E ^exp(^ 5 (x / ,x l ,f ;*))(!!) 

where S a is the classical action S of the trajectory 
a running from Xj to xy in time t, and Z?q = 
|det(5 2 5 5 /ax 4 9x / )| 1 /2 exp (-fp a ) is the Van Vleck de- 
terminant including the Maslov index. 

To characterize the decoherence process, we calculate 
the purity Tr (p|(i)), where the trace is performed over 
the particle degrees of freedom since ps is already a re- 
duced density matrix. If ps is initially pure, the purity 
starts from one and decays as time evolves due to heat 
bath influence. The semiclassical expression for the pu- 
rity is therefore given by the trace of the product of two 
ps{t) each one evolved by the semiclassical version of J 
calculated from (flOt when K is replaced by (fTTj) . That 
yiel ds the following expression which recalls the LE one 
[23l . [3(3] (we refer to Sec. V for a brief discussion about 
this issue) 



Tr(p 2 s W) 



(J- 

\2ttH 

x E 



dxfdyfdxidyidxIM J DfD S P[f]P[ S ] ps(x i ,y i ;0)p s (y^x^;0) 



J2 D & Dl,D n Dt, 
->y/.t). 

~(Sa(x/,Xi,f;*) - Si^y/j.J;*) + 5jj(y/,y^g;t) - SV(x/,X;,g;t) 



(12) 



a(x 4 ->x/,t), j)(y^y/,«), 

a'(yi^y/,t) jj'Cx^x/,*) 



x exp 



Due to the rapidly oscillatory phase factor contain- 
ing the action differences, most of the contributions will 
cancel out except for the semiclassically small action dif- 
ferences originating from pairs of trajectories which are 
close to each other in configuration space. We can thus 
use a linear approximation in order to relate the actions 
Sa, Sfj> along the trajectories a, fj' to the actions S a , S v > 
along nearby trajectories a, rj' connecting the midpoint 

= (xj + x^)/2 with X/. In the same way we relate 
Sa' > Sfj to the actions S a t, S n along the nearby trajecto- 
ries a', 77 connecting the midpoint r- = (y$ + y-)/2 with 
y/. The expansion will contain terms up to zero order 
in the D's and terms up to first order in the exponential. 



For Sa and S a > , the linearization yields [30j 

Sot ( x / ? x i if ) t) ~ 

r* 1 

S 0>a (xf,ri;t)- / dsq, a (s) ■ f (s) - -u p?, 
Jo 1 

Sf){yf,y'i,g;t) « 

/"* 1 

So,r,(yf,Ti,t) - J dsq v (s) -g(s) + -u'-p? 

(13) 

where S 0j0i is S [x] along the trajectory a, u = (Xj — xQ, 
u' = (y s ; — yQ and pf , are the initial momenta of 
trajectories a, 77. Analagous expressions are obtained for 
Sa' and Sfj' . As mentioned before, for the prefactors we 
have D a ps D a , Z)^ ps and analogously for D a >, D^i. 

In order to evaluate the sums over paths in (|12|) we 
consider only the diagonal contribution obtained by the 
pairings a — rj 1 and a 1 = r\. Besides that, we perform 
the gaussian averages over f and g as described in ©. 
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Using (JTJ) and ([8]) we obtain 

Tr(p|(i)) = J dxfdyfdndr'PsiruPDpJir'iiPi) 

x \ D «\ 2 \ D v\ 2 



(14) 



2k 



x exp 



/i 2 



J ds(q n {s) - q Q (s) 



where k = Am^ksT and 
^ (r ' P)= (2^) / du ^( r+ Y' r -? 



- ^;0] c-^ up 
(15) 



is the Wigner function of the initial ps- 



IV. DECOHERENCE RATES 

We consider two kinds of contributions in order to eval- 
uate (H3). First, if a and r\ are close enough to each other 
in phase space, it is possible to linearize the motion of 
one trajectory around the other to obtain 



qrj(s) - q«(s) 



I \ e As 

-t(p?-p?) V 

toA I 2 



(r' i -r i )--I-(p?-p?))^ 
mA / 2 



-As 



(16) 



where A is the Lyapunov exponent. On the other hand, 
if they are not close a and r\ can be first considered as 
free particle trajectories for time scales shorter than an 
average free flight time t Q . In this case, 

q„(V) - q«(«) « W - r<) + -(p? - pf ). (17) 

TO 

After t , we can only estimate (q. n (s) — q Q (s)) 2 by its av- 
erage value over the billiard area A. Assuming ergodicity, 
this average can be calculated as 



(.q„(»>-q„(*>r) A = / <W (q J ? ■ (18) 



It can be verified for simple geometries that (fT5|) yields 

H(s)-q«( s )f>i« A 

The prefactors \D a \ 2 and ID,,! 2 in (fl"4")) can be regarded 
as Jacobians when transforming the integrals over the 
final positions x/ and y/ into integrals over the initial 
momenta p^ = p» and p" = p- [H,l2l| leading to 

Tr(p|(t))= J dvdr'^dv'^ {^^pj W^) 



x exp 



2k r 

"fi2 J ds( 9 2 ( r i> r 'i)Pi)Pi; S ) 



(19) 



where the function Q 2 is taken either by the square of 
(|H>| or (fTT)) or simply by (|T£|) replacing p? by p^ and p" 
by p-. Thus, for initial states given by single gaussian 
wave packets as 



Ps I ' 
1 

(to 5 ) 



u u 

2' r -2 ; °) = 



exp 



(r-r ) 2 



I^ 2 



(20) 



Tr (p|(t)) can be calculated for the different situations 
mentioned above. 

Defining r = yt and using (JTTJ) , the result for times 
between and r Q = jt is 



Tr(p|(r)) = 



1 



16aiT + g a 2T 3 



1 + T T 



(21) 



where a\ — ksT/E, E = 



D = 



ik B T 
1717 



7f 
2auT- 



and 02 = -D/717 2 with 



Since our model is valid for times much 



shorter than I/7, r is certainly smaller than 1. The val- 
ues of a\ and 02 define relations between the free parame- 
ters of the problem. It is important to mentioned that we 
have not assumed at any point a highly localized initial 
state. Therefore, a 2 ~ 1 when a 2 is comparable to D/j 
and ai ^> 1 for localized wave packets. Eq. (|2"TT) shows 
that the purity starts from one, as it should be since the 
initial state (|20p is pure, and decreases as r approaches 

To- 

When t > t , there are two contributions: either the 
trajectories are close to each other in the sense mentioned 
before or they are not. In the latter case, since (fT5|) is 
mainly A, (JTIJ together with ([20]) yield 



(22) 



Tr (p 2 (r)) = exp 



1R keT < 

— l07T — T 

A V 



where A = is the mean level spacing for a billiard 
of area A. In the regime considered here, fc^T s> A and 
([22]) decays very fast. 

When the trajectories are correlated by the Lyapunov 
spreading, (fToT) and (gDJ yield 

Tr (p|(r)) = |l + h [(1 + 62) sinh (2Ar) + (1 - fe 2 )2Ar] 
cosh (2Ar) - - 1 J , (23) 

where b x = b 2 = (2E/HX) 2 , b 3 = 32 (£^£) 2 and 

A = A/7. For Ar > 1, only those terms with e 2Ar are 
relevant and (|2"3"]) can be written as 



Tr(p 2 s (r)) 



J2Kt\ 



(24) 



where (3 = h^^- + \. Eq. ^ shows that the classi- 
cally chaotic dynamics induces a decoherence rate inde- 
pendent of the heat bath parameters. 
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V. CONCLUSIONS AND DISCUSSIONS 

Summarizing, we have shown how to describe the de- 
coherence process of a classically chaotic system coupled 
to a heat bath within a path integral framework. Our 
approach leads to a decoherence rate given by the Lya- 
punov exponent. Although the results above were de- 
rived only for single gaussian wave packets with arbi- 
trary widths as initial states, the semiclassical approach 
presented also allows the treatment of their superposi- 
tions. This is an important difference compared to pre- 
vious works where onl y lo calized wave packets have been 
considered [l^, UUHil . The results for superpositions 
will be presented elsewhere and will be compared to the 
predictions of Refs.[l4|, EH which have claimed that the 
Lyapunov regime does not depend on the initial state. 
Concerning the description of the stable and unstable di- 
rections in (|16[) . that is not the most general one but 
it certainly captures the qualitative behavior one should 
find in a specific case. Different time scales of the classical 
chaotic dynamics were taken into account in the present 
calculation differently from those in Refs.[l4|, [H| where 
it is considered up to its linear regime only. 

It is possible to extend the present results beyond the 
high-temperature limit as long as the correlation time of 
thermal fluctuations is much shorter than the time scales 
of the system of interest dynamics Corrections to 
the diagonal approximation performed here could also 
be studied [27j for action differences of the order of h. 
As in the LE case [30| . they would take quantum effects 
into account systematically. Dissipation was neglected in 
the present calculation making the heat bath treatment 
comparable to previous ones using master equations [l4l . 
fl6j . However when the coupling to the heat bath is strong 
enough, dissipative terms cannot be neglected and it is 
still an open question whether the Lyapunov decoherence 
rate would be robust to that. 

Concerning the experimental observation of our theo- 
retical results, it was predicted recently in Ref.[3l| that 
Josephson junctions devices could also be used to ob- 
serve the Lyapunov regime in the time evolution of the 
fidelity. However, those devices are almost isolated from 
the heat bath in that context. We believe that the same 
regime and initial state preparation described there could 
be used to study the decay of interference fringes and to 
observe the Lyapunov decoherence rate as long as the 
heat bath temperature is increased. 

We briefly compare now the semiclassical calculation of 
the purity presented here and those of the LE fl9l [23l . [30ll 
At first glance, eq. tfTiZj) is identical to eq.(5.7) in Ref.|23j 



or to eq.(73) in Ref . (30| . They certainly have one thing in 
common: all of them are given in terms of four sums over 
chaotic trajectories which are considered unperturbed ei- 
ther by the heat bath influence (in the purity case) or by 
the extra potential (in the LE case). In other words, 
those equations result from the same perturbative ap- 
proach that allows us to use the well-known properties of 
those trajectories. One obtains then from those equations 
the leading order result after performing diagonal ap- 
proximation, which is another common point they have. 

The first subtle (though important) difference concerns 
the initial state. In Ref.[23[ for example, as most of the 
semiclassical analytical works on LE, a further approxi- 
mation is performed over eq.(5.7) since a highly localized 
gaussian wave packet is considered as initial state lead- 
ing to eq.(5.11). A more general calculation however was 
done recently in Ref. [30] allowing any kind of initial state 
and which we have followed closely here. The second dif- 
ference concerns the averages. Eq.(73) in Ref.[30J shows 
very clearly that LE is mostly calculated as an average of 
fidelity amplitude squared modulus. On the other hand, 
the purity is obtained from the trace of the reduced den- 
sity matrix squared, i.e. of the product of two averaged 
quantities. Hence the averages over the random potential 
that appear in (|12p are always performed independently. 
In LE case, the perturbation often depends on the tra- 
jectories which may be correlated or not. Thus averages 
cannot be performed always independently. Neverthe- 
less both purity and LE semiclassical expressions have 
two kinds of contributions arising from uncorrelated and 
correlated trajectories, the last one leading to the Lya- 
punov regime. The uncorrelated trajectories give rise to 
the Fermi golden rule regime in the LE case, which some- 
times dominates its decay. For the purity they lead to a 
very big decoherence rate (eq.([22])) that kills this contri- 
bution very fast and makes the Lyapunov rate always the 
dominant one (this result agrees with the numerical ones 
in Ref. [l6|). Finally, we have shown that the Lyapunov 
regime can be obtained from the exponents (eq. (|15[) and 
([19])) instead of the prefactors [H|,[23| avoiding short time 
problems. 
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